Astronomy & Astrophysics manuscript no. ChandraPMsastroph 


© ESO 2009 


February 5, 2009 





Proper motions of thermally emitting isolated neutron stars 

measured with Chandra 



C. Motch\ A. M. F. HaberP, A. Schwope"*, and V.E. Zavlin^ 



o 
o 

(N 
oo 

o 



> 

o 
o 



o 

O 



X 



' CNRS, Observatoire Astronomique, Universite de Strasbourg, 1 1 rue de TUniversite, F-67000, Strasbourg, France 

e-mail: motchOastro.u-strasbg . fr 
^ Instituto de Astronomia, Geoffsica e Ciencias Atmosfericas, Universidade de Sao Paulo, R. do Matao 1226, 05508-090 Sao Paulo, 

Brazil 

^ Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstrasse, D-85748 Garching, Germany 

Astrophysikalisches Institut Potsdam, An der Stemwarte 16, D-14482 Potsdam, Germany 
' Space Science Laboratory, Universities Space Research Association, NASA MSEC VP62, Huntsville, AL 35805, USA 

Received ; accepted 

ABSTRACT 

The remarkable astrometric capabilities of the Chandra Observatory offer the possibility to measure proper motions of X-ray sources 
with an unprecedented accuracy in this wavelength range. We recently completed a proper motion survey of three of the seven 
thermally emitting radio-quiet isolated neutron stars (INSs) discovered in the ROSAT all-sky survey. These INSs (RX J0420.0-5022, 
RX J0806.4-4123 and RXJ1308. 6-1-2127) either lack an optical counterpart or have one so faint that ground based or space born 
optical observations push the current possibilities of the instrumentation to the limit . Pairs of ACTS observations were acquired 3 to 5 
years apart to measure the displacement of the sources on the X-ray sky using as reference the background of extragalactic or remote 
Galactic X-ray sources. We derive 2cr upper limits of 123 mas yr"' and 86masyr"' on the proper motion of RX J0420.0-5022 and 
RX J0806.4-4123, respectively. RXJ1308. 6+2127 exhibits a very significant displacement (~ 9 cr) yielding fi = 220 + 25 masyr"', 
the second fastest measured among all ROSAT discovered INSs. The source is probably moving away rapidly from the Galactic 
plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of ~ 740 
((i/700 pc) km/s might be the largest of all ROSAT INSs and its corresponding spatial velocity stands among the fastest recorded for 
neutron stars. RXJ1308. 6-1-2127 is thus a middle-aged (age ~ 1 My) high velocity cooling neutron star. We investigate its possible 
origin in nearby OB associations or from a field OB star. In most cases, the flight time from birth place appears significantly shorter 
than the characteristic age derived from spin down rate. Overall, the distribution in transverse velocity of the ROSAT INSs is not 
statistically different from that of normal radio pulsars. 

Key words, stars: neutron. X-rays: general 



1. Introduction 

Observational studies of neutron stars started more than forty 
years ago with the discovery of radio pulsars. Owing to the 
large luminosity released by accretion, neutron stars in bina- 
ries were among the first sources detected by the pioneering 
X-ray instruments. The advent of sensitive imaging X-ray tele- 
scopes such as Einstein, ROSAT and now Chandra and XMM- 
Newton has revealed the complex pattern of high energy emis- 
sion from isolated neutron stars (INSs). More than forty radio 
pulsars have been detected in X-rays, with an X-ray luminos- 
ity ranging from 10"^ to 10"^ of their spin-down luminosity 
(Kargaltsev & Pavlov, 2008). Although radio surveys are still 
the main providers of new INSs, X-ray observations are con- 
tinuously revealing the existence of neutron stars undetected in 
the radio regime. These radio-quiet neutron stars are unlikely 
to constitute an homogeneous class. Among these are central 
compact objects (CCO) in SNRs. These soft X-ray sources are 
characterized by a steady thermal X-ray emission with a tem- 
perature in the range of 0.2 to 0.5 keV (see Pavlov et al. (2004) 
and de Luca (2008) for recent reviews). The ages of their as- 
sociated SNR are of a few kyr Two CCOs exhibit pulsations 
periods of a fraction of a second and their period derivative 
points at low surface magnetic fields of less than a few 10" G 
(Gotthelf & Halpern, 2008). Anomalous X-ray pulsars (AXPs) 



and soft-y repeaters (SGRs) form other classes of INSs revealed 
through their luminous X-ray emission. While SGRs appear as 
young as CCOs, their strong magnetic fields (B S 10''*G) drives 
a wealth of outstanding magnetospheric processes (see Woods 
& Thompson (2006) for a recent review). ROSAT was instru- 
mental in discovering the X-ray thermal emission of seven iso- 
lated neutron stars. Nicknamed "The Magnificent Seven", X-ray 
dim or X-ray thermally emitting isolated neutron stars (depend- 
ing on authors' preferences), we will simply refer to them as 
ROSAT discovered INSs. These INSs exhibit overall proper- 
ties that are at variance with those of most radio pulsars. Deep 
searches failed to detect any radio emission (Kondratiev et al., 
2008) although unconfirmed claims of detection at long wave- 
lengths exist (Malofeev et al., 2007). Their X-ray spectra are 
essentially thermal with kT in the range of 40 to 100 eV. Most 
ROSAT INSs also display shallow broad absorption features at 
energies around 300 eV which can be interpreted as proton cy- 
clotron lines or as atomic H or He transitions in high magnetic 
field conditions (> few 10'^^ G) (Ho & Lai, 2004). High proper 
motions in excess of 100 mas yr"' rule out accretion from the in- 
terstellar medium as a significant powering mechanism and thus 
leave no doubt that these INSs are young or middle-aged (age < 
10^ yr) cooling objects. Their X-ray spectra undergo little photo- 
electric absorption indicating distances of a few hundred parsec. 
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confirmed in two cases by HST parallaxes. With the exception 
of RX J 1605.3 +3249 for which no periodicity is detected yet, aU 
exhibit rather long spin periods in the range of 3.4 s to 11.4 s. 
When measured, the magnetic field implied by the spin down 
rate is roughly consistent with that inferred from the low energy 
absorption lines (see Haberl (2007) and Kaplan (2008) for recent 
reviews). 

Young neutron stars are expected to be strong X-ray sources, 
the high energy emission being either due to cooling or to non- 
thermal mechanisms. X-rays are thus an unbiased mean to detect 
nearby neutron stars. In this respect, it is striking that ROSAT 
discovered INSs appear to represent about half of all young (age 
< 3 X 10® yr) INSs known within ~ 1 kpc (Popov et al., 2003). 
ROSAT discovered INSs could thus be the tip of the iceberg of a 
formerly hidden large population of stellar remnants which may 
be, at least locally, as numerous as radio pulsars. 

Their properties lie at the crossroads between those of other 
groups of INSs. High magnetic fields and possible evidences 
for significant magnetospheric emission detected in the opti- 
calAJV bands may establish an evolutionary link between these 
objects and the group of magnetars (AXPs and SGRs) (Zane et 
al., 2008). Some long period and high B radio pulsars share the 
same locus as ROSAT discovered INSs in the PIP diagram, sug- 
gesting that the lack of radio emission could be due to the fact 
that the radio beam which narrows at long spin period does not 
sweep over the earth. Alternatively, the mechanism quenching 
radio emission could be similar to that transitively at work in 
rotating radio transients (RRATs) (McLaughlin et al., 2006). 

Proper motion measurements provide a number of important 
diagnostics on population properties. First, the relative velocities 
with respect to the interstellar medium put constraints on the 
contribution of accretion to the overall energy budget. Second, 
the reconstruction of backward trajectories allows to estimate 
the dynamical age of the neutron star assuming a birth in the 
Galactic plane and more precisely in one of the nearby OB as- 
sociations. The flight time can then be compared with the age 
derived from spin down, revealing a possible non standard brak- 
ing mechanism or pointing toward a long spin period at birth. 
Finally, the dynamical age also constrains the cooling mecha- 
nism (Yakovlev & Pe thick, 2004). 

Proper motions in the range of 108 to 333 mas yr ' have been 
measured for the three X-ray brightest ROSAT discovered INSs 
which have bright enough optical counterparts with niB S 27 
to allow measurements with large telescopes on the ground or 
with the HST (Kaplan et al., 2002b; Motch et al., 2003, 2005; 
Zane et al., 2006; Kaplan et al., 2007). In addition, the high as- 
trometric accuracy of the HST also allows to derive parallaxes 
for the nearest members (Kaplan et al., 2002b, 2007). However 
the remaining ROSAT INSs either have no known optical coun- 
terpart or have an optical candidate too faint (m^ > 28) to al- 
low repeated observations in a reasonable amount of time. X-ray 
imaging telescopes of the former generation had an insufficient 
spatial resolution to detect source displacements on the sky in 
a reasonable interval of time, even for high proper motion ob- 
jects such as INSs. The advent of high spatial resolution X-ray 
observatories has made possible the measurement with attract- 
ing accuracy of the proper motions of objects too faint to be 
detected at other wavelengths. Radio-quiet neutron stars are ob- 
viously prime targets for this kind of studies. We thus took ad- 
vantage of the unprecedented spatial resolution of the Chandra 
Observatory to start in 2002 first epoch measurements of three 
of flie X-ray faintest members, RX J0806.4-4123, RX J0420.0- 
-5022 and RXJ1308. 6+2127 with second epoch taking place 
three to five years later. Observing close to the maximum of the 



spectral energy distribution Chandra provides a fast and efficient 
way to measure accurate relative positions with respect to the 
background of active galactic nuclei. In this paper we report on 
the results of this observing campaign, extensively discuss the 
error budget on the proper motion and present the implications of 
our findings on our understanding of the properties of this class 
of neutron stars. Preliminary results were presented in Motch et 
al. (2007) and Motch et al. (2008). 

2. Target properties 

2.1. RXJ0806.4-4123 

The source was discovered by Haberl et al. (1998) in the ROSAT 
all-sky survey. XMM-Newton observations revealed pulsations 
at a period of 11 .4 s (Haberi & Zavlin, 2002; Haberi et al., 2004), 
the longest recorded so far for this group of objects. With a 
best fit blackbody temperature of yfeT = 96eV, RX J0806.4-4123 
lies among the hottest ROSAT discovered INSs. The spectral fit 
is significantly improved by adding a shallow absorption line 
centered around 460 eV. The interstellar absorption toward the 
source is Nh = 4xl0''cm~^ or l.lxlO^^'cm"^ depending on 
whether the low-energy broad absorption line is included or not 
in the spectral fit. RX J0806.4-4123 is the ROSAT INS located 
at the lowest Galactic latitude (b = -4.98°). Its modest photoelec- 
tric absorption suggests small distances of the order of 240 pc 
(Posselt et al., 2007). No optical counterpart brighter than a non 
constraining limit of ~ 24 is present in the error circle (Haberl 
et al., 1998). 

2.2. RXJ0420. 0-5022 

RXJ0420.0-5022 is flie X-ray faintest INS discovered in flie 
ROSAT all-sky survey (Haberi et al., 1999). With kT = 45 eV, 
RX J0420.0-5022 is the coolest of its group. It also displays a 
broad absorption line at E ~ 330 eV (Haberl et al., 2004). The 
source is located at relatively high galactic latitude (b - -44.39°) 
and the intersteUar absorption toward RXJ0420. 0-5022 is Nh 
= 1.0x10^° cm~^ or 2.0x10^" cm"^, again depending on the in- 
clusion or not of a broad line in the spectral model, and is on 
average about twice that of RX J0806.4-4123. XMM-Newton 
EPIC observations have revealed a pulsation period of 3.45 s, the 
shortest of all ROSAT discovered INSs. Haberi et al. (2004) also 
reported the possible existence of a = 26.6+0.3 mag object 
in the Chandra error circle. However, the presence of a relatively 
bright star close to the INSs adds a lot of scattered optical light 
at the source position and makes deep imaging difficult. The dis- 
tance estimate by Posselt et al. (2007) on the basis of interstellar 
absorption is of the order of 345 pc. 

2.3. RXJ1 308.6+21 27 

RXJ1308.6+2127 (alias RBS 1223) is the fifth X-ray brightest 
of the INSs discovered by ROSAT. It exhibits one of the highest 
temperatures of the group kT ~ 102eV (Schwope et al., 2007) 
and a pulse period of 10.3 1 s. A pair of possible harmonic proton 
cyclotron absorption lines are detected at ~ 0.23 and ~ 0.46 keV 
which would indicate a magnetic field of the order of 4x10'^ G 
(Schwope et al., 2007). However, the ratio of the line fluxes does 
not match that expected from cyclotron absorption and alterna- 
tive interpretations in terms of two distinct poles or atomic tran- 
sitions may also be considered. Kaplan & van Kerkwijk (2005b) 
could connect pulse phases observed by Chandra and XMM- 
Newton over a five year long time interval and derive a spin 
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down rate of P = 1.12x10 s/s. If due to magnetic dipole ra- 
diation P implies B = 3.4x10'^ G, a value roughly consistent 
with that estimated from the low energy absorption lines and 
the spin down age is ~ 1.5x10* yr. The optical counterpart of 
RXJ1308.6+2127islikelyafaintmv ~ 28.6 object Kaplan etal. 
(2002a) detected thanks to deep HST observations. The source 
is located at high galactic latitude (b = +83.08°). 

3. Observations 

Up to now, only few attempts have been made to measure proper 
motions in X-rays. The first example is the confirmation of the 
high proper motion of RX J1856.5-3754 (0.34 ± 0.12 masyr'^ 
Neuhauser (2001)) based on two ROSAT HRI observations ob- 
tained three years apart. A second case is the derivation of an up- 
per limit of ~ 170 mas yr"' on the motion of AXP 1E2259+586 
(Ogelman et al., 2005) using ROSAT, Chandra and XMM- 
Newton data. Another study bears on the central compact object 
in Puppis-A for which Hui et al. (2006) derive a proper motion 
of 107.5 + 34.5 mas yr"' by comparing two Chandra HRC-I ob- 
servations. Very recently, de Luca et al. (2008) and Kaplan et al. 
(2008) using Chandra AClS-1 observations reported 90% upper 
Umits on the proper motion of two magnetars, SGR 1900+14 
and IE 2259+586, of 54masyr~^ and 65 masyr"\ respectively. 

Our observing campaign started in Chandra Cycle 3 (2002) 
during which we obtained first epoch images for the three 
targets, RX J0806.4-4123, RX J0420.0-5022 and RXJ1308.6- 
+2127. Second epoch observations took place in 2005 and 2007 
(see Journal of Observations in Table 1). We selected the ACTS 
cameras because of their superior hard X-ray sensitivities com- 
pared to that of HRC-I. We also thought that the impossibility to 
filter X-ray events according to energy would possibly degrade 
the centering accuracy of X-ray sources with HRC-I. Finally, at 
the time we had to take the decision, the fact that HRC pixels 
are determined by the electronic readout was felt as being a risk 
for the long term stability of the astrometry. Our goal was to use 
the background of extragalactic sources, mainly AGN, to define 
an astrometric reference frame relative to which we could mea- 
sure the displacement of the isolated neutron star. In the case of 
the low Galactic latitude field of RX J0806.4-4123 {h — 4.98°) 
some Galactic active stars are also expected to contribute to the 
reference frame. However, being active stars (mostly younger 
than 1-2 Gyr ; Motch et al. (1997)), their proper motion is small. 
Stars younger than 1-2 Gyr still share most of the motion of their 
parent molecular clouds and thus have a small random velocity 
dispersion of less than ~ 20kms~'(Wielen, 1977). Being low 
luminosity X-ray sources (Lx^ lO^^ergs"'), the active coronae 
detected by Chandra are located at relatively short distances (~ 
1 kpc). Accordingly, the expected proper motion resulting from 
velocity dispersion and Galactic differential rotation is of the 
order of 22masyr-i at 300 pc for RX J0806.4-4123. Averaged 
over the whole ACIS-I field of view, the effect of young star dis- 
placements on relative astrometry is therefore negligible com- 
pared to the overall error budget. 

First MARX^ simulations carried out in 2001 indicated that 
an exposure time of 20 ks was sufficient to detect a suitable num- 
ber (> 10) of well localized background reference sources and 
obtain an astrometric error of ~ 0.07" on the relative positioning 
of the two epochs. Because of the enhanced soft X-ray sensitivity 
of ACIS-S compared to ACIS-I, the two X-ray brightest targets, 
RXJ0806.4-4123 and RX J1308.6+2127 would have suffered 
from pile-up in ACIS-S and would thus have likely produced 



Table 1. Journal of Observations 



Object 


Observation 


Camera 


Exposure 


Roll angle 




date 




time (ks) 


(deg) 


RXJ0806.4-4123 


2002-02-23 


ACIS-I 


17.7 


322.4 


RXJ0806.4-4123 


2005-02-18 


ACIS-I 


19.7 


325.6 


RXJ0420.0-5022 


2002-11-11 


ACIS-S 


19.4 


18.6 


RXJ0420.0-5022 


2005-11-07 


ACIS-S 


19.7 


23.6 


RXJ1308.6+2127 


2002-05-21 


ACIS-1 


19.5 


226.5 


RXJ1308.6+2127 


2007-05-12 


ACIS-1 


20.2 


226.7 



biased X-ray positions. Using ACIS-I for these two sources pro- 
vided the best compromise between the need to observe enough 
background sources and the requirement to have the best posi- 
tion for the INS target. However, for the fainter and much softer 
X-ray source RX J0420.0-5022, ACIS-I would have collected 
too few counts. Accordingly we chose to use ACIS-S at the ex- 
pense of a smaller field of view. The total source count rates were 
of 6x10-2, 6x10-2 and 1.1x10-^ ctss-' for RX J0806.4-4123, 
RX J0420.0-5022 and RX J1308.6+2127 respectively implying 
a pile-up fraction of less then 3% in the worst case. 

Our observing strategy was to enforce the same acquisition 
mode, roll angle and aimpoints at both epochs so as to minimize 
the impact of a different orientation on the sky or camera settings 
on the final accuracy of the relative astrometry. 

4. Data analysis 

4.1. Source detection 

Data reduction was performed in two steps. The 2002 and 2005 
observations of RX J0806.4-4123 and RX J0420.0-5022 were 
analysed using CIAO 3.3.0.1 while the 2002 and 2007 data of 
RXJ1308.6+2127 were reduced with CIAO 3.4 and CIAO 4.0. 
In all cases we used the most recently available calibration 
database. We corrected remaining attitude errors using the as- 
pect calculator available in the science threads^ at the Chandra 
X-ray Center (CXC). All level- 1 data were reprocessed to the 
level-2 stage using standard CIAO tasks. 

In order to remove the "gridded" appearance of the images 
and to avoid any aliasing effect, the standard Chandra data re- 
duction pipeline randomizes the X-ray event positions detected 
within a given pixel. This randomization process could slightly 
degrade the source centering accuracy. We thus removed the 
pixel randomization by reprocessing the raw data according to 
instructions available at CXC and analyzed these data in addi- 
tion to the usual randomized sets. 

Both the celldetect and wavdetect CIAO packages were used 
for source detection and characterisation. The task celldetect 
uses a shding cell of variable size to search for locations where 
the signal to noise ratio of a candidate source is larger than a 
given threshold. The position of the source is then derived from 
the center of gravity of the X-ray events found in the detect cell. 
In wavdetect, source detection is performed by convolving the 
X-ray image with wavelets of various spatial extents. The posi- 
tion of the source is eventually derived using the same algorithm 
(centroid) as for celldetect^ . 

Exposure maps were created in order to avoid spurious de- 
tections along the edges of the CCDs and obtain more accurate 
positions for those sources located close to these edges. These 



http://space.mit.edu/CXC/MARX/ 



^ http://a.sc. harvard.edu/ciao/threads/index.html 

^ http://cxc.harvard.edu/ciao/download/doc/detect_manual/ 
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maps were only used when running the celldetect task. Pixel 
scales of 1, 2 and 4 were used for wavdetect. These values are 
well suited to the detection of unresolved sources at moderate 
off-axis distances. We only considered sources detected in the 
central CCDs, chips 10 to 13 (16.9' by 16.9') and S2 and S3 (8.3' 
by 16.9') for ACIS-I and ACIS-S respectively. 

The final relative astrometric accuracy of the reference 
frames depends on the number of sources common to the two 
observations and on the quality of the determination of their po- 
sitions. We thus tried five different source detection thresholds, 
defined in celldetect as the significance of the source expressed 
in Gaussian sigma and in wavdetect as the probability of a false 
detection at any given pixel. These thresholds ranged from 1.5 
to 2.75 for celldetect and from 10"** up to a level of 5x10"^ 
for wavdetect. Testing several different thresholds allows to find 
the best compromise between the number of conmion sources 
and the mean positional errors, both rising with increasing de- 
tection sensitivity. The largest threshold value in wavdetect im- 
plies that about 40 false sources per CCD chip enter the source 
Ust. However, since we only consider sources detected at both 
epochs and located at a maximum distance of one to three arc- 
sec, the actual probability that the source common to both obser- 
vations is spurious is very low. Best positions are obtained in the 
energy range for which the signal to noise ratio of the reference 
sources is strongest. This range depends on the relative energy 
distribution of these sources and of the diffuse extragalactic and 
instrumental backgrounds. For the neutron star, measured in the 
same conditions as the reference sources, the choice of the en- 
ergy range is much less critical since the object is brighter than 
any of the background sources. We decided to test energy bands 
0.3-lOkeV, 0.3-5keV, 0.3-3keV, 0.5-5keV and 0.5-2keV. 

4.2. Matching reference frames 

The absolute astrometric frame derived independently for each 
Chandra observation is already very accurate since the 90% 
confidence radius is estimated to be of 0.6" close to aimpoint 
(Chandra X-ray Center). However, reaching the ultimate preci- 
sion on differential astrometry requires to correct for the remain- 
ing relative random attitude errors affecting the observations at 
the two epochs. We thus carried out a relative boresight correc- 
tion by allowing the reconstructed equatorial positions of the 
second epoch to be sUghtly shifted in right ascension and dec- 
lination and subjected to a small rotation around the aimpoint 
with respect to first epoch. We searched for translations of ± 1" 
in each direction with steps of 10 mas and for a rotation of up 
to + 0.1 deg with steps of 5x10"^ deg. The best offset and rota- 
tion angles were estimated using a maximum likelihood scheme, 
which for a bivariate normal distribution is equivalent to mini- 
mize the quantity 



Q 



n / 



Aai 



where Aq,, and A^, are the distances in right ascension and dec- 
lination between the positions of source / at epoch one and two 
after transformation. The errors on the corresponding distances 
are CTAai and ctasi and are computed from the individual source 
uncertainties in right ascension and declination returned by the 
centering algorithms. The accuracy with which we can link the 
astrometric reference frames of the two epochs is given by 



*^frame 



A — ^- — y(A2.+A2.) 
\ n(n-l) 



4.3. Systematic errors 

The ultimate accuracy with which Chandra can achieve relative 
astrometry is described in Chapter 5 of the Chandra Proposers' 
Observatory Guide. In this study, a sample of over 1100 sources 
with more than 50 counts detected in the framework of the 
Chandra Orion Ultra Deep Project (Getman et al., 2005) was 
boresight aligned and cross-correlated with the 2MASS cata- 
logue. After taking into account 2MASS positional errors, the 
distribution of the distances between ACIS-I and 2MASS en- 
tries indicated a systematic error with 90% confidence value 
of 0.15" on X-ray positions corresponding to a ~ 0.07" radial 
one sigma residual error. Therefore, although the celldetect and 
wavdetect algorithms can in principle center a relatively bright 
X-ray source on the CCD chip with a theoretical accuracy only 
restricted by signal to noise ratio, the ultimate positioning accu- 
racy is hmited by other systematic effects of unknown origin, but 
probably related to pixel to pixel sensitivity changes, intra pixel 
sensitivity variations, thermally induced misalignement, etc. We 
thus quadratically added a systematic error of o-gys = 70 mas to 
all ACIS source positions, including those of the bright target 
isolated neutron star. 

For our pairs of observations, o"frame is of the order of 
100 mas and thus dominates the error budget. The total error on 
the displacement of the central target is thus 



.2 

INS,\ 



-I- £7-? -I- 2 0-2, 
frame 



(3) 



5. Checking the actual accuracy achievable 

Only few attempts were done so far to measure proper motions 
in X-rays and the feasibility of such an endeavour is not of- 
ten discussed in the general literature nor in the technical doc- 
umentations. Note, however, that very recently Kaplan (2008) 
discussed to some extent the various errors entering astromet- 
ric AClS-1 observations and reached similar conclusions to ours. 
Before assessing the reality of the small displacement of an X- 
ray source on the sky one needs to envisage all potential sources 
of errors. This is best done by measuring in X-rays the motion of 
a source whose proper motion is accurately known from obser- 
vations at other wavelength, e.g. radio or optical. Unfortunately, 
to our knowledge, at the time of writing this paper, the Chandra 
database does not contain repeated ACIS observations of high 
proper motion isolated neutron stars with properties such that 
they can be used as test data. The only possibiUties left to test the 
error budget were then first, to perform extensive simulations of 
moving targets and second, to use repeated observations of non 
moving objects to test systematic errors. We examine these two 
options in the following subsections. 



(1) 5.1. Simulations 



(2) 



The goal of these simulations was twofold. First, to estimate the 
accuracy with which a small displacement of the INS could be 
detected using the ACIS-I and ACIS-S configurations, allow- 
ing for different aimpoint positions and systematic offsets of the 
reference frames. Second, to find the most efficient centering 
method and the best suited parameters, such as energy range, 
detection threshold and randomization state. 

We used MARX 4.2.1 to create fake event lists mimicking the 
2002 and 2005 observations of the fields of RX J0806.4-4123 
and RX J0420. 0-5022. The simulated observations had the same 
duration as the original ones and in addition, the same starting 
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times. This last constraint allows to account for the slow degra- 
dation of the quantum efficiency of the CCDs, especially at ener- 
gies below 1 keV. For each camera, the simulated sources were 
those common (within 3") to the actual 2002 and 2005 obser- 
vations as detected in the 0.5-5.0 keV band and with a threshold 
of 10-'^ for RX J0806.4-4123 and in the 0.5-2.0 keV band with 
a threshold of lO'^ for RX J0420.0-5022. A total of 26 and 12 
sources were created for ACIS-I and ACIS-S simulations respec- 
tively. Each simulated source had the same position on the sky 
and photon flux as the actual one. Since most of the background 
reference sources are too faint to constrain their energy distribu- 
tion, we assumed that they were all of extragalactic origin and 
that their spectra were described by an absorbed power law of 
spectral index F = 1.7, undergoing a hydrogen column density 
of A^H = 5.21 X lO^i cm-2 (for the field of RXJ0806.4-4123) 
and A^H = 1.07 x lO^" cm'^ (for the field of RX J0420.0-5022). 
The total Galactic column densities were derived from the far 
infrared emission maps of Schlegel et al. (1998) available at the 
NASA/IPAC Extragalactic Database and applying the Predehl & 
Schmitt (1995) relation between optical and X-ray absorptions. 
The energy distribution of the neutron stars were taken as black- 
bodies with kT - 96 eV and kT - 44 eV absorbed by Nh - 
4x10''^ cm-2 and Nh =1x10^'^ cm^^ for RX J0806.4-4123 and 
RX J0420.0-5022 respectively. Input count spectra were gener- 
ated using XSPEC 12.2.1. We created a diffuse background using 
the observed deep blank fields provided by the Chandra calibra- 
tion database, by re-projecting them to the tangent plane of the 
observations and for each simulation by randomly picking back- 
ground events so as to reproduce the level of background noise 
of the actual observation in the 0.5-5 .0 keV energy band. Finally, 
sources and background were merged into a single event file. 

We verified that simulated sources displayed a count his- 
togram consistent with the expected Poisson distribution. The 
input flux of each simulated source was modified in order to ac- 
count for the exposure map, in particular for those sources close 
to the edge of a CCD. We also checked that on average, the count 
rates of the simulated sources were identical to those of the ob- 
served ones. 

Table 2 lists the parameters of the 13 difl'erent settings used 
for the first and second epoch simulations. We applied three dif- 
ferent offsets to the position of the central soft X-ray source rep- 
resenting the neutron star in the event file of the second epoch. 
A number of parameters were changed as well. The position 
of the aimpoint was moved by a fraction of a CCD pixel in 
order to produce slightly different distributions of the source 
PSF on the pixel grid. Note, however, that the slow Lissajous 
shaped dithering applied to Chandra pointing during observa- 
tions should smear out PSF effects to a large extent. A system- 
atic offset of 0.29 arcsec (±0.15 "and ±0.25"in RA and DEC 
respectively) was applied to the 2005 coordinate system of run 
R2 to test our capability to correct for small shifts in the astro- 
metric frames. Additionally, the randomization in pixel - which 
is part of the standard processing in the Chandra pipeline - was 
disabled in an extra sample of simulations in order to investi- 
gate the impact of this effect on the final positional accuracy. 
All these artificial observations were eventually processed by the 
same pipeline as the real data using the two CIAO detection al- 
gorithms, celldetect and wavdetect. 

None of the two algorithms, celldetect and wavdetect could 
recover aU of the reference sources entered in the simulations, 
as a consequence of the background and photon counting noise 
introduced. At the maximum sensitivity used here, celldetect (cr 
= 1.5) finds less sources than wavdetect (threshold of 5x10"^). 
Averaging over all simulations and energy bands used, celldetect 



Table 3. Matching error (in arcsec) and detection method for ACIS-S 
and ACIS-I simulated observations. 



Field 




Camera 


Method 


(y frame 


rms 


RXJ0420.0- 


-5022 


ACIS-S 


celldetect 


0.164 


0.032 


RXJ0420.0- 


-5022 


ACIS-S 


wavdetect 


0.184 


0.038 


RXJ0806.4- 


-4123 


ACIS-I 


celldetect 


0.134 


0.022 


RXJ0806.4- 


-4123 


ACIS-I 


wciydetect 


0.128 


0.019 



finds on average 6.2 and 8.1 sources common to both epochs 
for the ACIS-S and ACIS-I simulations respectively while for 
wavdetect the figures are 7.0 and 14.8. 

We evaluated the relative merits of the two detection and po- 
sitioning methods using as criterion the statistical error derived 
on relative frame matching between the two epochs. Table 3 lists 
the distribution of these frame errors for ACIS-I and ACIS-S, 
averaged over energy bands and detection thresholds for which 
the number of common sources was comprised between 7 to 
10 for RXJ0420.0-5022 and 12 to 16 for RX J0806.4-4123. 
None of the two algorithms proves a clear-cut advantage, their 
performances remain statistically consistent with each other. 
Removing pixel randomization while keeping all other param- 
eters identical does sUghtly improve source centering accuracy. 
The effect remains however relatively small, at the level of a few 
percents. 



Not unexpectedly, the energy band considered for source de- 
tection also affects to some extent the accuracy with which X- 
ray sources are positioned. Wide bands allow more sources to 
be detected, but introduce a higher background. For our simu- 
lations, the wider (0.3-lOkeV) and narrower (0.5-2keV) bands 
tend to yield slightly worse errors on positions. The small off- 
sets (± 0.15" and ± 0.25") applied to the coordinate system of 
the second epoch of run R2 are also properly recovered. 

We finally checked the accuracy with which the small mo- 
tion of the simulated neutron star could be sized after frame 
matching. We show in Fig. 1 the measured displacement of the 
neutron star compared to the one entered in the simulation. We 
applied here to the reductions a selection similar, although not 
exactly identical, to that applied to the real data. Namely, we 
impose the presence of at least 7 reference sources for ACIS- 
S and 12 for ACIS-I and we discard reductions yielding frame 
errors larger than the mean of all reductions. Randomized and 
de-randomized simulations are included. Error bars show the 
dispersion of the results from the selected reductions and sim- 
ulations. Shifts as small as 0.34" are properly recovered with 
an accuracy of 0.138" and 0.075" for the ACIS-S and ACIS-I 
configurations respectively using the wavdetect algorithm. 

5.2. Repeated fields 

Our simulations gave us high confidence in the possibility to 
measure small displacements of unresolved X-ray sources with 
ACIS. The MARX ray tracing suite reUes on the characteristics 
of the telescope and cameras as measured before launch in the 
laboratory and later evolved to account for actual flight perfor- 
mances. However, a number of unadvertised effects might intro- 
duce additional errors on the final relative astrometry recover- 
able in real observations. We thus searched the Chandra archive 
for repeated observations of extragalactic targets fulfilling simi- 
lar conditions as ours in order to confirm that the relative posi- 
tioning errors were indeed consistent with those estimated from 
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Table 2. Simulation parameters. A RA and A Dec are the displacements applied to the aimpoint or to position of the neutron star in the second 
observation with respect to those in the first one. 



Run 


Epoch 


Aimpoint 


ARA 


A Dec 


Displacement 


ARA 


A Dec 






offset 


aimpoint 


aimpoint 


target 


target 


target 






(pixel) 


(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


Rl 


2002 


Nominal 


0.000 


0.000 










2005 


Nominal 


0.000 


0.000 


1.00 


-0.707 


-0.707 


R2* 


2002 


Nominal 


0.000 


0.000 










2005 


Nominal 


0.000 


0.000 


1.00 


-0.707 


-0.707 


R3 


2002 


0.25 


-0.123 


0.000 










2005 


0.25 


-0.123 


0.000 


1.00 


-0.707 


-0.707 


R4 


2002 


0.50 


0.000 


-0.246 










2005 


0.50 


0.000 


-0.246 


1.00 


-0.707 


-0.707 


R5 


2002 


0.67 


0.232 


0.232 










2005 


0.67 


0.232 


0.232 


1.00 


-0.707 


-0.707 


R6 


2002 


Nominal 


0.000 


0.000 










2005 


Nominal 


0.000 


0.000 


0.75 


-0.375 


-0.650 


R7 


2002 


0.25 


-0.123 


0.000 










2005 


0.25 


-0.123 


0.000 


0.75 


-0.375 


-0.650 


R8 


2002 


0.50 


0.000 


-0.246 










2005 


0.50 


0.000 


-0.246 


0.75 


-0.375 


-0.650 


R9 


2002 


0.67 


0.232 


0.232 
















n 9^9 


U. / J 






RIO 


2002 


Nominal 


0.000 


0.000 










2005 


Nominal 


0.000 


0.000 


0.34 


-1-0.232 


-1-0.253 


Rll 


2002 


0.25 


-0.123 


0.000 










2005 


0.25 


-0.123 


0.000 


0.34 


-1-0.232 


-1-0.253 


R12 


2002 


0.50 


0.000 


-0.246 










2005 


0.50 


0.000 


-0.246 


0.34 


4-0.232 


-1-0.253 


R13 


2002 


0.67 


0.232 


0.232 










2005 


0.67 


0.232 


0.232 


0.34 


-fO.232 


-hO.253 



* For this run, all source positions in the second epoch were offset by a fraction of an arcsecond (see text) 

simulations. Constraints were the existence of a rather bright un- 
saturated unresolved X-ray source with enough relatively faint 
reference sources around in pairs of ACTS images obtained with 
as far as possible similar roll angles and exposure times. Not 
surprisingly, very few observations match these conditions since 
only programmes aimed at long or short term monitoring are 
good candidates, the rest of the repetitions usually implying dif- 
ferent modes or different instruments. We eventually found 3 
pairs of ACIS-I and 3 pairs of ACIS-S observations suitable 
for the tests. We list in Table 4 the main properties of these 
fields. These pairs of observations were processed with exactly 
the same pipeline as that used for our astrometric programme. 
In four instances, the source whose position was monitored was 
the original AGN target. For the two "Field" observations we 
had to select a relatively bright source (with a total number of 
counts larger than ~ 40 in the shortest exposure) and located as 
close as possible to the aimpoint. As for our INSs observations, 
we estimated the displacement of the target source using the five 
combinations of detection thresholds and energy bands provid- 
ing the lowest crf^^e- We only discuss here the results of the 
wavdetect algorithm since celldetect yields quite similar frame 
errors and displacements. Removing pixel randomization does 
not necessarily improve frame matching. We list in Table 5 the 
measured displacements of the target source and their statistical 
significance. None of the extragalactic targets shows a signifi- 
cant motion on the sky since in all cases its significance remains 
below ~lcr. The largest shift is seen on the LALA Cetus field 
target which is also the X-ray faintest of our test sources. In 
all other cases, the displacements of the bright source measured 
after matching the astrometric frames is smaller than 100 mas. 
These results agree quite well with those of our simulations and 



Table 5. Measured displacement A (mas) of the central source in re- 
peated fields and its statistical significance. 



Field A cr A cr 



randomized no randomization 



LALA Cetus Field 


173 


0.92 


96 


0.59 


PKS 0312-7707 


100 


0.88 


112 


1.02 


FIELD- 142549-1-353248 


54 


0.38 


84 


0.58 


3C 33 


62 


0.50 


91 


0.70 


3C 325 


66 


0.57 


43 


0.40 


MG J0414-F0534 


112 


1.03 


48 


0.49 



give us further confidence in the astrometric capabilities of the 
Chandra Observatory. Interestingly, small shifts in the coordi- 
nates of the aimpoint and large changes in the roU angle do not 
seem to impact much the accuracy of the relative astrometry. 



6. Results 

Some reference sources showed rather large positional differ- 
ences between the two epochs. Although these position offsets 
were not statistically significant considering the associated er- 
rors, we investigated how the inclusion of these sources in the 
transformation from one epoch to the other was impacting the 
value of the frame errors. Several test runs performed on the 
actual data showed that best frame errors were obtained when 
only sources located at a maximum distance of 1" were included 
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Fig. 1. Mean displacements in a and d, averaged over the selected simulations (see text), compared to their input values (diamonds). Error bars 
show the rms of the measurements values. Black and red squares denote celldetect and wavdetect respectively. 




-1 -0.5 0.5 ] -1 -0.5 0.5 



Displacement in Right Ascension (arcsec) Displacement in Right Ascension (arcsec) 



Table 4. Repeated Chandra observations used as test fields. 



Target 


Observation 


Camera 


Mode 


Exposure 


Roll angle 


Target Used 




date 






time (ks) 


(deg) 




LALA Cetus FIELD 


2003-06-13 


ACIS-I 


VFAINT 


160.1 


121.5 


Source at 


LALA Cetus FIELD 


2003-06-15 


ACIS-I 


VFAINT 


14.2 


121.6 


02:04:46.5s +05:04:04.4 


PKS03 12-770 


1999-09-08 


ACIS-I 


VFAINT 


12.8 


63.8 


PKS 0312-770 


PKS03 12-770 


1999-09-08 


ACIS-I 


VFAINT 


12.6 


63.8 




FIELD- 142549+353248 


2002-04-16 


ACIS-I 


VFAINT 


120.1 


163.4 


Source at 


FIELD- 142549+353248 


2002-06-09 


ACIS-I 


VFAINT 


58.5 


219.9 


14:25:47.1 +35:39:54.8 


3C 33 


2005-11-08 


ACIS-S 


FAINT 


19.9 


281.45 


3C 33 


3C 33 


2005-11-12 


ACIS-S 


FAINT 


19.9 


281.45 




3C 325 


2005-04-14 


ACIS-S 


VFAINT 


29.6 


145.22 


3C 325 


3C 325 


2005-04-17 


ACIS-S 


VFAINT 


28.7 


145.22 




MG J04I4+0534 


2001-11-09 


ACIS-S 


FAINT 


28.4 


57.35 


MG J0414+0534 


MG J04I4+0534 


2002-01-08 


ACIS-S 


VFAINT 


96.7 


295.24 





in the transformation. Sources showing larger positional differ- 
ences were in almost all cases located at large distances to the 
optical axis. In other instances, one of the two sources or both 
sources had a very small detection significance suggesting a spu- 
rious detection. We also tried both the celldetect and wavdetect 
source detection and centering tasks. For ACIS-I the error on the 
transformation delivered by celldetect was 10 to 70% larger than 
that resulting from wavdetect positions. These results are at vari- 
ance with the outcome of our simulated runs and of the repeated 
fields which showed that the two methods provided very com- 
parable accuracies. In the case of ACIS-S (RX J0420.0-5022), 
however, the two algorithms yielded similar results. We also note 
that the value of the displacement of the central source does not 
vary significantly with the task used. Therefore, in the following, 
we will only consider results obtained with the wavdetect task. 

As for simulated data and repeated fields, the accuracy with 
which the two astrometric frames can be matched together de- 
pends on the randomization being applied or not, on the energy 
band used for source detection and on the value of the detection 



threshold. We show in Fig. 2 how the frame eiTor varies with 
the number of sources common to both epochs. This number de- 
pends on source detection threshold and on energy band. 

Not unexpectedly, the randomization process tends to 
slightly deteriorate source positioning, although in the case of 
RXJ0806. 4-4123, the situation is inverted. Better frame errors 
are also obtained on the randomized data of some of the repeated 
fields studied above. In a similar manner, the energy band pro- 
viding the best frame eiTor slightly varies with target. In general, 
the restricted 0.3-5 keV, 0.3-3.0 keV or 0.5-5 keV bands yield the 
smallest errors. This behaviour is consistent with the increasing 
ACIS background above 5 keV and lower number of photons in 
the narrower 0.5-2 keV band and is similar to that found in sim- 
ulations. 

Fig. 3 shows the average number of sources detected in each 
observation and the number of sources common within 1" to 
both observations as a function of detection threshold. Data are 
presented here for the 0.3-5 keV energy band. The threshold used 
in wavdetect represents the probability of a false detection at any 
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Fig. 2. Accuracy of the astrometric matching (frame error in arcsec) as a function of the number of common sources plotted for different energy 
bands. 
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Frame error 



' 

0.05 
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given pixel of the image. Using wavelet scales of sizes 1, 2 and 4 
pixels, over the central four ACIS-I CCD and central two ACIS- 
S chips, the expected number of false detections is of the order 
of 40 in ACIS-I at threshold = 10"^. As expected, the num- 
ber of detections dramatically rises with increasing thresholds 
reflecting the enhancement of the number of false sources. For 



RX J0806.4-4123 and RX J1308.6-F2127 (ACIS-I), the number 
of sources detected in each field increases by 24 when the thresh- 
old rises from 10"^ to 10"^. In contrast, the number of sources 
common to the two observations, which therefore have a high 
probability to be real increases much more slowly with thresh- 
old. 
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Fig. 3. Mean number of detected sources (solid Une) and number of 
sources common to the two observations within l"(dashed line) in the 
0.3-5 keV band. 




1 lO-'' 1 

Wavdetecl threshold 



Table 6. Measured frame offsets and rotation angle between two 
epochs. All distances are given in mas. 



Target 


RA offset 


DEC offset 


Rotation 




(mas) 


(mas) 


angle (deg) 


RXJ0806.4-4123 


156 


-398 


-0.021 


RXJ0420.0-5022 


248 


132 


-0.034 


RXJ1308.6+2127 


82 


266 


-0.010 



A significant number of high likelihood sources are not re- 
covered from one observation to the next. At very low thresholds 
of 10"^ or 10 for which the number of spurious detections 
should be completely negligible, only one third of all sources 
are detected at both epochs. Selecting the brightest sources only, 
still half of them are not recovered from one observation to the 
other. This strongly suggests that a large fraction of extragalactic 
sources, mostly AGN, are significantly variable on a time scale 
of several years. 

Overall, the density of field sources around the two INSs ob- 
served with ACIS-I are quite comparable and considering the 
very similar effective exposure times can be directly estimated 
from Fig. 3. The slightly higher number of detections around 
RXJ0806. 4-4123 could be due to a contribution by additional 
stellar coronae which owing to the very low Galactic latitude 
of the target are expected to be well detectable. The number of 
ACIS-S detections is also very close to half that in ACIS-I in 
accord with the ratio of the field of view used for the two config- 
urations. 

In order to be less sensitive to a particular configuration of 
reference sources and obtain more conservative errors, we aver- 
aged results obtained for the five (among a total of 25) combi- 
nations of detection thresholds and energy bands which yielded 
the lowest frame matching errors. Best matches of the astromet- 
ric reference frames were obtained with the parameters listed in 
Table 6. 



We list in Table 7 the displacements and their associ- 
ated errors measured between the two epochs and corrected 



for the values listed in Table 6. The positions of RXJ0806.4- 
-4123 and RXJ0420.0-5022 in 2002 and 2005 are consistent 
within one Gaussian sigma. The two sigma upper limits on the 
proper motions of RX J0806.4-4123 and RX J0420.0-5022 are 
86masyr ' and 123 mas yr ' respectively. 



In contrast, a very significant displacement is detected 
for RXJ1308.6-I-2127, with a statistical significance of ~ 9cr 
and corresponds to a total proper motion of fi - 220 ± 
25 masyr The magnitude and direction of the position move 
of RXJ1308.6-I-2127 with respect to the background of extra- 
galactic sources does not depend significantly on the randomiza- 
tion step being applied or not in the processing of the individual 
photons, nor does it depend on the number of reference sources 
considered. We also checked that the celldetect algorithm was 
yielding similar results. Expressed in J2000 equatorial coordi- 
nates, the proper motion vector is; 

fia cos(5) - -207 + 20 mas yr"' ; /ij = 84 ± 20 mas yr~^ 

and corresponding proper motions in Galactic coordinates 
are: fii cos(l) - -122 mas yr"' and /i/, = 183 masyr"'. The source 
is therefore moving toward larger Galactic latitudes as would be 
naively expected for a neutron star bom in the Galactic plane. 
We note, however, that for such a high Galactic latitude ob- 
ject, the velocity perpendicular to the Galactic plane is domi- 
nated by the unknown radial component. RX J1308.6+2127 dis- 
plays the second largest proper motion among the seven ROSAT 
discovered INSs, just after the X-ray brightest ROSAT INS 
RXJ1856. 5-3754. Several distance estimates are proposed in 
the literature. Kaplan et al. (2002b) derive a distance of ~ 670 pc 
for RXJ1308.6-I-2127 using the preliminary HST parallax of 
RX J 1856.5 -3754 (d~ 140 pc) and making the simplifying as- 
sumption that the optical emission arises from neutron star at- 
mospheres with similar surface compositions and areas. The re- 
vised distance of 160 pc for RX J1856.5-3754 (van Kerkwijk & 
Kaplan, 2007) increases the estimated distance of RX J 1308.6- 
-1-2127 to ~ 760 pc. Under these assumptions, the distance pre- 
dicted for RXJ0720.4-3 125, J ~ 340 pc, is in fair agreement 
with that derived from the HST parallax, d - 360^™ pc (Kaplan 
et al., 2007). Based on detailed X-ray light curve and SED mod- 
elling Schwope el al. (2005) derived distances in the range of 
76 to 380 pc. The short distance resulted for the extreme model 
which assumed small, isothermal hot spots on the cool stellar 
surface which remained undetected in X-rays. While this model 
explained the SED and the X-ray variability, it was regarded 
unlikely due to the unknown heating mechanism of the spots. 
The more realistic temperature model based on calculations by 
Geppert et al. (2004) resulted in a comparative representation of 
the data but a much larger distance. This number is subject to 
uncertainties given the uncertain temperature profile and the un- 
known stellar parameters. The distance of 380 pc was based on 
Rns - 12 km, Tgr = 4, hence = 13.8 km. The growth of the in- 
tegrated interstellar absorption with increasing distance has been 
used by Posselt et al. (2007) to constrain the distance of several 
ROSAT discovered INSs. This method applied to INSs having 
astrometric parallaxes yields consistent values. However, the lo- 
cal absorption models which rests on the measurement of the Na 
I line equivalent width in hot or fast rotating stars has a too sparse 
coverage at high Galactic latitude to be used effectively and has 
to be replaced by a less accurate analytical description. In ad- 
dition, the exact value of the photoelectric absorption derived 
from X-ray modelling somewhat depends on the number and 
position of the low energy shallow absorption lines. Fitting the 
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Table 7. Measured displacements. All values given in mas. 



Target 


Target error 


Frame Error 


Total error 


Measured 
displacement 


Significance 
(cr) 


RXJ0806.4-4123 


12 


81 


128 


117 


0.91 


RXJ0420.0-5022 


26 


151 


184 


209 


1.13 


RXJ1308.6+2127 


8.5 


73 


123 


1096 


8.9 



XMM-Newton spectra of RXJ1308. 6+2127 with two low en- 
ergy lines gives Nh ~ 1.2-1. 8xlO^°cm"^ (Schwope et al., 2007). 
This Nh value suggests a distance of at least ~ 500 pc (Posselt 
et al., 2007). In this respect, the very short distance of 76 pc pro- 
posed by Schwope et al. (2005) seems incompatible with the 
observed photoelectric absorption. The best guess distance esti- 
mate of RX J1308.6+2127 is thus probably in the range of -400 
to ~800pc. We note that the proper motions published in Motch 
et al. (2007) and Motch et al. (2008) derived from a less exten- 
sive analysis are fully consistent with those reported here. 

7. The possible age and birth place of 
RXJ1 308.6+21 27 

The accuracy with which Chandra is able to measure proper mo- 
tion vectors does not exactly compare with what can be achieved 
using ground based telescopes or the HST. Nevertheless, the 
present determination turns out to be quite useful to estimate the 
possible origin of the neutron star and its likely dynamical age. 

We computed the past trajectories of the neutron star dur- 
ing the last 4 Myr assuming a + 2cr range in fia and fig, present 
distances from 50 to 1200pc and radial velocities in the inter- 
val of ± 700 kms '. We corrected the observed proper motion 
vector for the displacement of the Sun with respect to the Local 
Standard of Rest. Considering the large velocity, young age and 
proximity of the neutron star and the relatively short time during 
which we follow it, we neglected effects due to Galactic grav- 
itational potential and differential Galactic rotation. For most 
trajectories, these effects should be of ^ 1% (Bienayme et al., 
2006). 

As already quoted above, a displacement toward higher 
Galactic latitude does not necessarily imply that RX J 1308. 6- 
-1-2127 escapes from the Galactic plane since atb = 83.082°, the 
unknown radial velocity dominates the component perpendicu- 
lar to the Galactic plane. As a matter of fact, any radial velocity 
toward the Earth in excess of a few tens of kms ' would im- 
ply that RX J1308.6H-2127 is born at large distances from the 
Galactic plane and the vast majority of the trajectories escaping 
from the plane require receding radial velocities. 

We first investigated whether some backwards trajectories 
would cross any of the known nearby OB associations. For 
that purpose we compiled a merged list from the catalogues in 
Humphreys (1978), Ruprecht et al. (1981) and de Zeeuw et al. 
(1999), keeping for each association the most accurate informa- 
tion available, mainly the boundaries and mean distances deter- 
mined by Hipparcos for the nearest groups. 

The number of past flight paths crossing any given OB as- 
sociation is low, less than ~ 2%. This reflects the small vol- 
ume of the targeted associations compared to that sampled in 
total, which is necessarily large, owing to the substantial uncer- 
tainties on the proper motion vector, current distance and more 
importantly on the unknown radial velocity. In addition, only 
few nearby OB associations are counted in the I ~ 0°-20° range 
of Galactic longitude where RX J1308.6-H2127 seems to origi- 



Table 8. Possible birth places of RX J1308.6-I-2127 listed with the 
present distance of the respective OB association and corresponding 
neutron star flight times, current distance and radial velocity ranges. 



OB assoc 


Distance 


Age range 


do 


Vr 




(pc) 


(Myr) 


(pc) 


(km/s) 


UpperSco 


145 


0.55+0.25 


260+50 


-430+180 


SctOB2 A 


510 


0.90+0.15 


535+53 


-550+95 


SctOB2 B 


1170 


1.38+0.26 


800±90 


-520±110 



nate from. The times at which backwards trajectories cross any 
given group of massive stars show a well defined peaked dis- 
tribution with a small tail extending to the oldest ages consid- 
ered in our study (4 Myr) caused by few very low radial veloc- 
ity trajectories. Hereinafter, current distances (do) and radial ve- 
locities (Vr) will refer to their mean over all trajectories com- 
prised within one full width at half maximum of the peak of the 
age distribution. Our analysis shows that only three OB associ- 
ations (see Table 8) could have given birfli to RX J1308.6+2127 
if we assume that the current distance to the neutron star is less 
than ~ 900 pc. Because of its proximity (d = 145 pc), the Upper 
Scorpius OB association (/ = -17°,0°; = 10°,30°) exhibits a sig- 
nificant motion which has to be taken into account when com- 
puting intersection times. The Upper Scorpius moving group is 
part of the general Sco OB2 association. Its past position was 
computed using group proper motion and radial velocity derived 
from Hipparcos and listed in de Zeeuw et al. (1999). According 
to Reichen et al. (1990), the Scutum OB2 association is likely to 
be the superposition on the line of sight of two disctinct groups, 
one located at 510pc and another at 1170pc. Intersections with 
more remote (d > 1.5 kpc) OB associations such as Sgr OBI, Set 
0B3, Sgr OB7, Ser OB2, Sgr 0B6, Ser OB 1 , Sgr OB4, Sgr 0B5 
all imply unreaUstic present distances to the source larger than 
do = 900 pc with average crossing times of the order of 1 .5 to 2.2 
Myr. 



However, a significant number of early type stars stars are 
found far from identified OB associations. The relative fraction 
of these field stars is not very well constrained and depends 
somehow on the accuracy of the photometric distances used and 
on the definition of the chosen boundaries for the associations. 
According to Garmany (1994), the frequency depends on spec- 
tral type, the fraction of early type field stars varying from ~ 
40% for stars hotter than B2 to 25% for O stars only. The major- 
ity of field objects are runaway stars identified from their high 
velocities (Blaauw, 1961; Stone, 1991) while the rest cannot be 
related to any group and might therefore be true field stars born 
outside clusters. Some of the runaway stars can convincingly be 
associated to the birth of radio pulsars (Hoogerwerf et al., 2001) 
and were thus kicked off their parent associations in the super- 
nova explosion in a massive binary. Alternatively, the existence 
of many binary runaway stars (Gies & Bolton, 1986) brings sup- 
port to the dynamical ejection scenario in which early gravita- 
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Fig. 4. Histograms of the crossing times of the Galactic plane for vari- 
ous assumed present distances to the source. 
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tional interaction with one or more cluster stars is responsible 
for the jump in velocity (Poveda et al., 1967). 

In a second step, we thus relaxed the condition of a birth 
place close to a catalogued OB association and only imposed an 
origin in the supernova explosion of a field star assumed to be 
located within the Galactic plane. 

Fig. 4 shows the histograms of the ages at which the 
backwards trajectories intercept the Galactic plane for different 
ranges of assumed present distances to the source. In a similar 
manner as for trajectories impacting OB associations, for a given 
range of current distances, the age histogram of all trajectories 
crossing the Galactic plane displays a clear peak with a reduced 
tail extending to longer times. Although ages older than 2 Myr 
appear possible for some parameters, the most probable ages 
peak at values increasing from 0.4+0.2 Myr for do - 100-400 pc 
to 1.3+0.3 Myr for d/o = 800-1000 pc (the interval of ages quoted 
here corresponds to the full width at half maximum of the peak 
of the histogram). For the most probable "intermediate" present 
distance range of 400 to 800 pc, the age distribution peaks at 
0.96 + 0.30 Myr. Galactic longitudes at birth span a sector from 
approximately / = -4° to / = +24°. Old birth dates correspond 
to large present distances and more remote birth places ranging 
from 130pc (do = 100-400pc) to 1600 pc (do = 800-1200pc). 

Note that the random position of the progenitor star with re- 
spect to the Galactic plane smears to a small extent the travel 
times computed above. The scale height of O to B5 stars be- 



longing to the local Galactic disk is ~ 40-60 pc (see Elias et al. 
(2006) and references therein). Taking this effect into account 
adds a scatter of ~ 10^ yr to the ages shown in Fig. 4. 

The main parameter determining the age of RXJ1308.6- 
+2127 is its present distance which is therefore the most im- 
portant source of uncertainty. A best guess do of 400 to 800 pc 
leaves as most likely birth place candidate the OB association Set 
OB2 A and an age of 0.90+0.15 Myr. If the early type star from 
which RXJ1308. 6+2127 originates was a field object in the 
Galactic plane, then comparable ages of 0.96 + 0.30 Myr appear 
as most probable. However, the much closer distances proposed 
by Schwope et al. (2005) allow for a birth in the nearby Upper 
Scorpius part of the Sco OB2 complex 0.55±0.25 Myr ago. 
RXJ1308. 6+2127 could thus be like RX J1856.5-3754, also 
born in the Upper Scorpius (Walter, 2001), and as RX J0720.4- 
-3125 likely bom in the Trumpler 10 or Vela OB2 associations 
(Motch et al., 2003; Kaplan et al., 2007), another product of the 
Gould Belt. 

8. Discussion 

The proper motion of RXJ1308. 6+2127 is the second largest of 
all ROSAT discovered INSs (see Table 9). However, owing to 
its plausible larger distance than other high proper motion mem- 
bers of the group, in particular that of the closest ROSAT INS 
RX Jl 856.5-3754, its tangential velocity Vt is probably consid- 
erable, in the range of 420 to 840km s"^ Backwards trajectories 
crossing either a catalogued OB association or more generally 
the Galactic Plane also imply rather large radial receding veloci- 
ties of the order of 550 km s \ a value weakly dependent on the 
assumed present-day distance. The 3-D velocity of RX J 1308. 6- 
-1-2 1 27 might thus be of the order of 600 to 1 000 km s~ ' . Such an 
enormous space velocity is nevertheless not unusual for neutron 
star standards. The radio pulsar B1508-I-55 is a good example of 
a hyper-fast neutron star (Vx = 1083;^^^"' km s"') where both as- 
trometric proper motion and trigonometric parallaxes were mea- 
sured (Chatterjee et al., 2005). The Guitar Nebula and its asso- 
ciated radio pulsar B2224+65 (Cordes et al., 1993) is another 
example of ultra high velocities (Vt = 1640 ((i/1900pc) km s"'). 
The high velocity of RXJ1308. 6+2127 precludes accretion of 
matter from the tenuous interstellar medium present at the large 
distance above the Galactic plane (do ^ 400 pc) where the source 
is currently located (see e.g. Treves et al. (2000)). This INS is 
thus without doubt another young cooling neutron star. 

Our Chandra observations of RX J0806.4-4123 and RX- 
J0420. 0-5022 yield significant constraints on the tangential ve- 
locities of these two other sources. At a distance of 345 pc de- 
rived from an analysis of the distribution of interstellar absorp- 
tion on the line of sight (Posselt et al., 2007), the transverse 
velocity of RX J0420.0-5022 already ranks among the lowest 
for the group. Moreover, the upper hmit of 86 mas yr"^ obtained 
on the yearly displacement of RX J0806.4-4123 is well be- 
low those displayed by the four sources for which a measure- 
ment exists. The probable small distance of the source (d ~ 
235 pc; Posselt et al. (2007)) implies comparatively low trans- 
verse velocities of less than 96kms"'. This brings the possibil- 
ity that RX J0806.4-4123, being also the source located at low- 
est Galactic latitude among all ROSAT INSs, could power part 
of its X-ray emission through accretion of matter from the in- 
terstellar medium. This heating mechanism seems however un- 
Ukely since its X-ray spectrum, surface magnetic field and spin 
period do not difl^er much from those of other high velocity mem- 
bers of the group. Besides, the high magnetic field B ~ 6.1 - 8.6 
xlO^^ G inferred by Haberi et al. (2004) from the energy of the 
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Table 9. Proper motions, distances and transverse velocities of the seven 
ROSAT discovered INSs. Upper limits are at the 2 a level 



Name 




Distance 


Vt 


Ref 




(mas yr"') 


(pc) 


(km s"') 




RXJ1856.5-3754 


333±1 


3611^ 


254 


1,2 


RXJ0720.4-3125 


107.8±1.2 


184 


3,4 


RXJ1605.3+3249 


155.0+3.1 


-390 


286 


5,6,7 


RXJ1308.6+2127 


220.3+25 


400-800 


417-835 


8,7,9 


RXJ2143.0+0654 




-430 




7 


RXJ0806.4-4123 


<86 


-235 


<96 


7,9 


RXJ0420.0-5022 


<123 


-345 


<200 


7,9 



References: 1) Kaplan et al. (2002b); 2) van Kerkwijk & Kaplan (2007) 
3) Motch et al. (2003); 4) Kaplan et al. (2007); 5) Motch et al. (2005) 
6) Zane et al. (2006); 7) Posselt et al. (2007); 8) Schwope et al. (2005) 
9) this work 
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Fig. 5. Transverse speeds of ROSAT discovered INSs compared to those 
of non recycled radio pulsars. Data taken from Hobbs et al. (2005). 
The black thick line and the thin magenta line show the histograms of 
the entire population and of the young (age < 3 xlO*yr) respectively. 
Measured values and upper limits for the ROSAT INSs are shown on 
the top. 



shallow absorption lines present in the soft X-ray spectrum and 
the rotation period of 1 1 .37 s would make accretion from a typ- 
ical low to medium density interstellar medium impossible (see 
e.g. Motch (2001)). In addition, the lack of long term variability 
of the X-ray flux is at variance with what would be expected for 
an accreting neutron star. 

The present work adds three new measurements or upper 
limits and doubles the size of the sample of ROSAT INSs for 
which information on proper motion exists, leaving the seventh 
member, RXJ2143. 0-1-0654 for further studies. It is now possi- 
ble to better compare the distribution in transverse velocities of 
the ROSAT discovered INSs with that of radio pulsars. 

We plot on Fig. 5 the distribution of the transverse velocities 
of non recycled radio pulsars, total and young (age < 3xlO''yr) 
populations as given in Hobbs et al. (2005). The positions of the 
six ROSAT INSs are also shown for comparison. ROSAT INSs 
exhibit 2-D velocity distributions consistent with those of radio 
pulsars. A simple Kolmogorov-Smirnov test gives a probability 



larger flian 50% fliat flie distribution of the ROSAT INS tangen- 
tial velocities is identical to that of the young (age < 3xlO''yr) 
radio pulsars. RXJ1308. 6+2127 fills in the large velocity tail 
of the radio pulsar distribution thus strengthening the similar- 
ity between the two populations. Within the limits enforced by 
small number statistics, it seems that the different birth condi- 
tions leading to the large magnetic fields and long spin periods 
of the ROSAT INSs are not related to a specific kick mechanism. 
This extends the apparent lack of dependency of radio pulsar 
properties on kick velocity (see e.g. Wang et al. (2006)) toward 
larger B and spin period. 

RXJ1308.6H-2127 is the fifth brightest INS in the ROSAT 
PSPC band and might be the first to be born outside of the Gould 
Belt, either in one of the Set OB associations or from a field 
OB star This INS could thus mark the approximate flux bound- 
ary beyond which the production of other more distant OB as- 
sociations or field stars becomes noticeable. We note that the 
0.1-2.4keV PSPC count rate of 0.29ctss"' is roughly consis- 
tent with the population model of isolated cooling neutron star 
by Posselt et al. (2008) which predicts that at PSPC count rates 
below ~0.1ctss ' INSs located beyond the Gould Belt should 
dominate. 

Assuming that the distance to RXJ1308. 6+2127 is indeed 
smaller than 800 pc yields ages of less than 0.80 Myr (for a birth 
in Upper Sco), 1.05 Myr (for a birth in Set OB2 A) or 1.26 Myr 
for a birth from a field early type star. In contrast, the charac- 
teristic age derived from spin down tsd -1.5 Myr is somewhat 
longer. Similar and even greater disagreements between dynami- 
cal and spin-down ages are found for the two other ROSAT INSs 
for which P was measured. In the case of RX J0720.4-3125, the 
flight time from the most likely parent OB associations is in the 
range of 0.5 to 1 Myr for a characteristic age of 1 .9 Myr (Kaplan 
et al., 2007). Recently, van Kerkwijk & Kaplan (2008) reported 
the determination of the spin down of RXJ1856. 5-3754 with 
again a huge discrepancy between the dynamical age of 0.4 Myr, 
consistent with the estimated cooling age of ~0.5Myr and the 
much longer spin down age of 4 Myr. 

Several mechanisms have been put forward to explain this 
discrepancy (see Kaplan & van Kerkwijk (2005a) and Kaplan 
& van Kerkwijk (2005b)). Among these are a long birth period 
of the order of several seconds and an early additional braking 
due to interaction with a fallback disc. The last mechanism en- 
visaged is magnetic field decay. Originally proposed by (Heyl & 
Kulkarni, 1998) it has been recently investigated in greater de- 
tails by Aguilera et al. (2008). First, field decay provides through 
Joule heating an additional source of energy which accounts for 
a slower decrease of temperature with time. Second, the charac- 
teristic times derived from P then strongly overestimate the true 
age of the neutron star Both effects converge to reconcile the rel- 
atively high polar temperatures observed in RX J0720.4-3125 
and RX J 1308. 6+2 127 with dynamical ages of up to 1 Myr while 
explaining the discrepancy between dynamical and spin-down 
ages. In addition, the strongly mass dependent cooling efficiency 
makes the age temperature relation even more complex. If both 
mass and dynamical ages were known with enough accuracy, it 
could become possible, in principle, to estimate the strength of 
the magnetic field at birth and in the case of ROSAT discovered 
INSs find out whether these objects have evolved from a magne- 
tar population (age -10* yr) or if middle-aged, were born with 
lower fields in the range of lO^^-IO'^^G (Aguilera et al., 2008). 
The unusually large optical excess detected in RX J2143. 0+0654 
(Zane et al., 2008) together with the large magnetic field evi- 
denced from the soft X-ray absorption lines (B ~ lO'^^G; Zane 
et al. 2005) hint at a close resemblance and perhaps evolution- 
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ary link between this ROSAT INS and magnetars. The magnetic 
field of RX J1308.6+2127 (B ~ 4x101^ G; Schwope et al. 2005) 
ranks among the lowest of the ROSAT INSs (Haberl, 2007). 
Decay from a typical magnetar field (B > few lO'^^G) would 
require an interval of time of the order of Tohm ~ 1 Myr. A large 
present distance to the source could thus bring some support to a 
magnetar origin for RX J1308.6+2127. However, the estimated 
relative birth rate of ROSAT INSs and magnetars seems to ex- 
clude that the entire population of ROSAT ESTSs evolves from 
former magnetars (Popov et al., 2006) 

9. Conclusions 

Chandra Observatory ACIS data obtained at a time interval of 
three and five years allow us to constrain the proper motion of 
RXJ0806.4-4123 and RXJ0420.0-5022 with two sigma up- 
per limits of yu < 86 mas yr~' and 123 mas yr"' respectively. We 
measure the very significant displacement of RXJ1308. 6+2127 
implying ^ - 220.3 + 25 mas yr ' . The present work doubles the 
number of radio-quiet thermally emitting neutron stars for which 
information on proper motion exists. Extensive simulations car- 
ried out with MARX and analysis of fields observed several times 
by Chandra confirm the feasibility of measuring proper motions 
at X-ray wavelengths with an accuracy approaching that of opti- 
cal ground based observations. 

The likely transverse velocity of RX J1308.6+2127 is large, 
of the order of 400 to 800 km s"' and could well be the largest of 
the six ROSAT discovered INSs for which a measurement of the 
proper motion and an estimate of the distance both exist. The re- 
ceding radial velocities required for a birth close to the Galactic 
plane are of the order of 500 km s"' suggesting space velocities 
as high as 1,000 kms"'. RX J1308. 6+2127 is thus another ex- 
ample of a young cooling radio-quiet nearby neutron star. Our 
upper limit on the displacement on the sky of RX J0806.4-4123 
hints at a slow relative velocity with respect to the interstellar 
medium. We argue, however, that accretion is unlikely to occur 
and does not significantly contribute to the X-ray emission. 

Within the limits imposed by the small number statistics, 
the distribution in transverse velocity of the ROSAT INSs does 
not appear statistically different from that of the radio pulsars. 
This indicates that the mechanism leading to the strong mag- 
netic fields and long periods typical of these objects is at first 
order independent from that controlUng the kick velocity. 

The assumed present distance to RX J1308.6+2127 deter- 
mines its possible travel time from nearby OB associations and 
from the Galactic plane in general. The analysis of all possible 
backwards trajectories suggest that for a present distance com- 
prised between 400 and 800 pc, the most Ukely birth OB asso- 
ciation is the Scutum OB2 A group located at 510pc whereas a 
birth in the nearby Upper Scorpius part of the Sco 0B2 associa- 
tion would be possible if the present distance to the source were 
much shorter {d ~ 260 pc). We also find that travel times from 
the Galactic plane or from one of the candidate OB associations 
depend sensitively on the assumed current distance to the source 
and range from less than ~ 1 .2 Myr for do ~ 400 - 800 pc to less 
than 0.6 Myr for shorter distances. Although the still rather large 
uncertainties on the proper motion vector, current distance and 
the unknown radial velocity allow for longer flight times in some 
cases, the majority of the trajectories imply a kinematic age sig- 
nificantly smaller than the spin-down time of 1.5 Myr. A sim- 
ilar discrepancy occurs for RX J0720.4-3125 and RX J1856.5- 
-3754. 
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